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SUMMARY 


The  reciprocal  equations  for  relating  the  incompressible  circula¬ 
tory  indicial  lift  to  the  lift  due  to  harmonic  oscillations  have  been 
modified  to  include  the  noncirculatory  lift  associated  with  apparent-  - 
mass  effects.  Although  the  apparent-mass  effects  are  impulsive  in 
nature  in  incompressible  flow,  the  lift  due  to  apparent-mass  effects  in 
compressible  flow  is  a  time-dependent  function.  The  corresponding 
reciprocal  equations  for  the  total  compressible  lift  are  given.  By  use 
of  the  reciprocal  equations  for  compressible  flow,  the  indicial  lift  and 
moment  functions  due  to  an  airfoil’s  experiencing  a  sudden  acquisition 
of  vertical  velocity  are  determined  numerically  for  Mach  number  0. 7. 

Lack  of  sufficient  flutter  coefficients  prevents  the  calculation  of 
these  functions  at  other  Mach  numbers. 

Although  the  indicial  lift  and  moment  functions  due  to  penetration 
of  a  sharp-edge  gust  may  be  obtained  from  the  oscillatory  tab  or  aileron 
coefficients  by  a  similar  analysis,  sufficient  coefficients  are  not 
available  at  the  present.  However,  an  approximate  method  is  shown  for 
determining  a  portion  of  this  unsteady-lift  function. 

When  a  comparison  is  made  of  the  indicial  lift  functions  at  Mach 
numbers  0.0  and  0. 7>  it  is  noted  that  the  growth  of  lift  to  the  steady 
state  appears  to  be  less  rapid  for  the  compressible  case  than, for  the 
incompressible  case.  Consequently,  the  calculation  of  the  gust  load 
factor  at  high  subsonic  Mach  numbers  utilizing  the  two-dimensional 
incompressible  indicial  lift  functions  and  an  over-all  correction  for 
compressibility  such  as  the  Prandtl-Glauert  factor  might  be  conservative. 


INTRODUCTION 


A  knowledge  of  transient  flows  is  important  in  many  aeronautical 
problems.  In  the  study  of  transient  flows,  two  types  of  airfoil  motions 
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have  had  special  significance  -  a  harmonically  oscillating  airfoil  and 
an  airfoil  experiencing  a  sudden  change  in  angle  of  attack.  The  lift 
function  for  an  airfoil  experiencing  a  sudden  change  in  angle  of  attack 
and  the  lift  function  associated  with  the  growth  of  lift  on  an  airfoil 
due  to  penetration  of  a  sharp-edge  gust  are  commonly  referred  to  as 
indicial  lift  functions.  The  present  paper  is  concerned  with  the  use 
of  reciprocal  relations  for  determining  for  compressible  flow  the 
indicial  lift  functions  directly  from  the  lift  data  that  are  available 
for  the  airfoil  oscillating  harmonically. 

The  indicial  lift  functions  have  been  derived  for  two-dimensional 
incompressible  flow;  the  function  for  a  sudden  change  in  angle  of 
attack  was  derived  by  Wagner,  see  reference  1,  while  the  penetration 
function  was  derived  by  Kussner,  reference  2.  An  account  of  the  rela¬ 
tions  that  exist  between  these  indicial  lift  functions  and  the  lift 
coefficients  for  a  two-dimensional  oscillating  airfoil  is  given  by 
Garrick  in  reference  3* 

In  a  recent  paper  by  Lomax,  Heaslet,  and  Sluder  (reference  4),  a 
method  for  determining  the  indicial  lift  and  moment  function  is  given. 
While  the  beginning  portion  of  the  indicial  functions  can  be  calculated 
readily  and  the  final  value  is  considered  to  be  the  steady-state  lift 
value  given  by  the  Prandtl-Glauert  factor,  the  intermediate  or  transi¬ 
tion  values  of  the  indicial  lift  functions  are  difficult  to  obtain  by 
this  method  and  consequently  numerical  results  are  given  only  for  a  Mach 
number  of  0.8  in  reference  4. 

The  situation  in  regard  to  the  subsonic  compressible-flow  coef¬ 
ficients  for  an  oscillating  airfoil  is  much  better.  Possio  (refer¬ 
ence  5)  has  formulated  the  problem  in  a  linearized  form  for  determining 
the  lift  and  moment  for  the  oscillatory  case.  Dietze,  Schade,  and 
Frazer  and  Skan  (references  6,  7j  and  8)  present  the  lift  coefficients 
of  a  harmonically  oscillating  airfoil  at  various  Mach  numbers  up  to  0.8 
and  for  various  values  of  reduced  frequency  up  to  about  2.5*  It  was 
felt  desirable  to  see  whether  these  available  flutter  data  could  be  used 
in  conjunction  with  the  reciprocal  relations  to  obtain  the  complete 
indicial  lift  functions  for  compressible  flow  since,  for  dynamic-load 
studies,  knowledge  of  the  indicial  lift  functions  is  needed  over  a 
larger  range  of  chord  lengths  than  that  given  in  reference  4  for  Mach 
numbers  other  than  0.8.  This  paper  discusses  the  use  of  the  reciprocal 
relations  for  the  case  of  compressible  flow  and  presents  an  evaluation  of 
the  indicial  functions  for  M  =  0.7  for  an  airfoil  suddenly  acquiring  a 
vertical  velocity,  for  which  case  it  was  found  that  sufficient  oscillatory 
lift  data  were  available. 

The  indicial  lift  function  due  to  penetration  of  a  sharp-edge  gust 
could  also  be  determined  by  superposition  provided  the  oscillatory  lift 
coefficients  were  known  for  a  wide  range  of  flap  to  chord  ratios  as  well 
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as  for  a  range  of  reduced  frequencies.  Unfortunately,  sufficient  data 
are  not  available  at  the  present  time;  however,  in  order  to  provide 
some  insight  as  to  the  growth  of  lift  for  the  gust  case,  an  alternate 
method  is  shown  for  approximating  most  of  this  function. 


SYMBOLS 


s,  s^,  a  distance  traveled,  half-chords 


CD 

V 

c 

k 

L(s) 

M(s) 

P 

h 

.  •• 

h,  h 

M 

8(b) 


angular  frequency 

forward  velocity  of  airfoil 

chord 

reduced- frequency  parameter 
lift  per  unit  length  of  span 

moment  per  unit  length  of  span  about  quarter- chord  point 
density 

amplitude  of  vertical  displacement,  half- chords 
first  and  second  derivatives  of  h  with  respect  to  s 
Mach  number 

impulse  function,'  6(s)  =  00  for  s  =  0,  8(s)  =  0  for  s  f  0, 

and  /  b(s)  ds  =  1 


k-^s) 


U 

k2(s) 


indicial  lift  function  for  an  airfoil  experiencing  a'  sudden 
acquisition  in  vertical  velocity  as  used  in  equation 

L  =  -7tpcV2hk1(s) 
velocity  of  sharp-edge  gust 

,  indicial  lift  function  due  to  penetration  of  a  sharp-edge 
gust  as  used  in  equation  L  =  -jrpcVUkgCs) 


k 
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m-^(s)  indicial  moment  function  for  an  airfoil  experiencing  a 

sudden  acquisition  of  vertical  velocity  where  moment  is 
taken  about  quarter-chord  position 

xCp(s)  center- of-pres sure  location  in  percent  chord  from  leading 
edge  for  an  airfoil  experiencing  a  sudden  acquisition  of 
vertical  velocity 

C(k)  Theodorsen*s  circulatory  lift  function  (reference  9)> 

F(k)  +  iG(k),  as  used  in  equation  L  =  -jtpcV2heiks  Q-ikC(k7| 

Cc(k)  complex  compressible-flow  oscillatory  lift  coefficient  which 

includes  both  circulatory-  and  noncirculatory-lift  com¬ 
ponents,  Cc(k)  =  Fc(k)  +  iGc(k),  as  used  in  equation 

L  =  -*pcV2eiksh£ikCc(k)] 

M(k)  in-phase  component  of  the  complex  moment  corresponding  to 

the  Fc(k)  lift  coefficient 

f(k)  =  Fc(k)  -  Fc(oo) 

m(k)  =  M(k)  -  M(°°) 

Zp,  Z2  in-phase  and  out-of-phase  lift  coefficients  associated  with 
translation  of  the  airfoil  as  used  in  equation 

L  =  jtpcV2elks  |^ZX  +  1Z2) 

Mp,  M2  in-phase  and  out-of-phase  moment  coefficient  about  quarter- 
chord  position  associated  with  translation  of  the  airfoil 

as  used  in  equation  M  =  jtpc2V2eiks  +  iM2) 

A,  B,  C  constants 
Matrix  notation: 


rectangular  matrix 


* 


*r 

W 


* 

* 
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ANALYSIS  AND  EESULTS 

* 

Superposition  Integrals  for  the  Incompressible 
and  Compressible  Cases 


The  reciprocal  equations  which  are  essentially  a  special  form  of 
the  superposition  integral  have  been  derived  for  the  incompressible 
circulatory  lift  by  Garrick  (reference  3).  This  section  first  discusses  - 
an  extension  of  Garrick’s  work  by  incorporating  the  potential  lift  due 
to  apparent-mass  effects  into  the  reciprocal  equations  for  the  incom¬ 
pressible  case.  The  reciprocal  equations  for  the  compressible  case  are 
then  indicated  and  the  evaluation  of  the  indicial  lift  and  moment  for 
an  airfoil  experiencing  a  sudden  acquisition  of  vertical  velocity  is 
determined  for  the  case  of  M  =  0.7. 

Incompressible  case.-  In  the  incompressible  case,  the  lift  in 
unsteady  motion  has  been  separated  into  two  parts  -  a  noncirculatory 
lift  which  is  commonly  referred  to  as  an  apparent-mass  effect  and  a 
lift  due  to  the  circulation  about  the  airfoil.  The  reciprocal  equa¬ 
tions  for  incompressible  flow  are  expressed  only  in  terms  of  Theodorsen’s 
circulatory- lift  function  C(k)  and  Wagner’s  indicial- lift  function  k-jjs) 
and  can  be  expressed  in  two  separate  forms  as  follows: 


ki(s)=i/0  —  fn  ks  (s>0)  (la) 

and 

ki<s>  -it§/0“G(k)  r kSd*  (b>o)  <ib> 


where  F(k)  and  G(k)  are,  respectively,  the in-phase  and  out-of-phase 
components  of  circulatory  lift  on  a  harmonically  oscillating  airfoil  defined 
by  the  equation 


"  L  =  -itpcV^heiks[likC(k)] 

% 


C(k)  =  F(k)  +  iG(k) 


where 
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<T 


The  k-^s)  function  is  defined  such  that  the  lift  on  an  airfoil  experi¬ 
encing  a  sudden  acquisition  of  vertical  velocity  may  he  expressed  hy  the 
following  equation: 


L  =  -jtpcV^^s)  (2) 


Two  relations  that  can  he  obtained  from  equation  (la)  are  (l)  that 
the  value  of  F(k)  at  k  =  0  is  equal  to  the  asymptotic  value 
of  k-]_(s)  as  s  tends  to  infinity,  these  values  corresponding  to  the 
steady-state  lift,  and  (2)  that  the  initial  value  of  k-^(s)  is  equal 

in  magnitude  to  the  asymptotic  value  of  F(k)  as  k  tends  to  infinity. 

This  second  relation  is  not  readily  evident  from  equation  (la);  con¬ 
sequently,  a  mathematical  proof  of  this  second  relation  is  shown  in 
appendix  A.  As  a  result  of  these  relationships,  the  initial  value  and 
the  asymptote  of  the  k^(s)  function  can  he  determined  from  the  given 

F(k)  function.  A  knowledge  of  these  end  points  is  valuable  in  a 
numerical  solution  of  the  reciprocal  equations. 

As  noted  previously,  equations  (la)  and  (lb)  express  only  the 
circulatory  part  of  the  lift  and,  therefore,  do  not  account  for  the  * 

apparent-mass  effects.  An  expression  for  the  apparent- mass  effects  for 
the  indicial  case,  however,  may  he  obtained  as  follows. 

* 

The  lift  due  to  apparent-mass  effects  alone  for  any  arbitrary 
vertical  motion  of  the  airfoil  -when  written  in  terms  of  the  nondi- 
mensional  displacement  h  is  given  in  reference  9  as 


L  ‘  -”p  T I  S(s> 


(3) 


When  this  equation  is  applied  to  the  case  of  an  airfoil  experiencing  a 

sudden  vertical  velocity,  h(s)  becomes  impulsive  in  character.  The 
magnitude  of  this  impulse  may  be  defined  by 


h(s)  =  h&(s) 


where 


» 

h 


is  the  instantaneous-  vertical  velocity  acquired  by  the  airfoil 


* 
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and  the  function  &(s)  is  defined  such  that  8(s)  =  ~  for  s  =  0, 
S(s)  =  0  for  s  f  0,  and 


&(s)ds 


When  h5(s)  is  substituted  for  the  impulsive  sinking  acceleration  h(s) 
in  equation  (3),  the  resulting  expression  for  the  lift  due  to  apparent 
mass  alone  becomes 


^impulse  ~  -itpcV^h  -  (1).) 

If  the  impulsive  part  of  the  unsteady-lift  functions  is  designated  by 
kl(s) impulse  an<i  tlie  form  of>  equation  (2)  is  retained,  then  the 
expression  for  the  impulsive  lift  may  be  written 


w 


•'■'impulse 


-jtpcV  hk^( s) impulse 


(5) 


The  value  of  k-^s)  due  to  apparent-mass  effects  in  the  incompressible 

case  may  then  be  determined  by  a  comparison  of  equations  (4)  and  (5)  to 
give 


kl^s^impulse  2  (6) 

Addition  of  the  k-^(s)  function  due  to  circulation,  equations  (l),  and 
the  ki(s)injpulBe  function,  equation  (6),  then  gives  the  expressions 
for  the  lift  for  the  indicial  case  as  follows: 

kl(s)total  =  |  kS  &  +  i  6(s)  (s  £  0)  (7a) 
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and 


9 


^l( 8 ) total 


(s  =  0)  (7b) 


Compressible  case.-  While  in  the  incompressible  case  it  has  been 
found  convenient  and  natural  to  consider  the  circulatory  and  non- 
circulatory  lifts  separately  (this  separation  is  especially  desirable 
because  of  the  impulsive  nature  of  the  noncirculatory  part),  these 
circumstances  do  not  exist  in  the  compressible  case.  The  perturbation 
velocities  are  finite  in  the  compressible  case  (since  speed  of  sound  is 
considered  finite)  and,  consequently,  the  effects  of  the  pressure  dis¬ 
turbances  in  the  vhole  flow  field  are  not  felt  instantaneously  as  for 
incompressible  flow.  Consequently,  the  lift  due  to  apparent-mass 
effects  would  not  be  impulsive  but  would  remain  finite  and  time  dependent. 
Separation  of  the  lifts  is  not  necessary  in  the  compressible  case  and, 
moreover,  is  not  convenient  since  the  available  results  for  the  com¬ 
pressible  case  are  presented  numerically  only  for  the  total  lift.  The 
reciprocal  equations  for  the  total  compressible  lift  are,  therefore, 
much  simpler  than  the  equations  for  the  incompressible  case  since  the 
impulse  function  may  be  disregarded.  w 

For  compressible  flow,  the  case  is  exactly  analogous  to  the  treat¬ 
ment  of  indicial  admittance  given  by  Bush  in  reference  10.  In  this  y 

reference,  equations  are  given  which  relate  the  indicial  admittance  to 
the  in-phase  and  out-of-phase  responses  of  a  system  subjected  to  an 
oscillatory  forcing  function.  In  notation  similar  to  that  used  for 
incompressible  flow,  these  expressions  may  be  written  as  follows 


kl( s) total 


Fc(k)  sin  ks 
k 


dk 


(s  >0)  (8a) 


and 


kl(s) total 


Fc(0) 


2  r°°  Gc(k)  cos  ks 

n  Jq  k 


dk 


(s  >  0)  (8b) 


where  Fc(k)  and  Gc(k)  are,  respectively,  the  in-phase  and  out-of-  * 

phase  lift  components  on  an  oscillating  airfoil  and  are  defined  by  the 

equation  v 
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L  =  -itpcV2helks  -ikCc(k)j 


(9) 


where 


Cc(k)  =  Fc(k)  +  iGc(k) 


Equations  (8a)  and  (8h)  are  applicable  for  the  compressible  case  pro¬ 
viding  the  function  Cc(k)  is  continuous  and  finite  in  the  interval  0 

to  oa.  From  the  physical  conditions  for  the  compressible  case,  it  can 
reasonably  be  assumed  that  these  functions  adhere  to  those  conditions. 

It  is  of  interest  to  note  that,  when  these  equations  are  applied  to 
incompressible  flow,  they  apply  only  to  the  circulatory  part  of  the  lift 
as  given  by  equations  (la)  and  (lb).  In  this  case,  the  factor  Fc(0) 

becomes,  by  definition,  simply  unity.  For  compressible  flow,  the  fac¬ 
tor  Fc(0)  acquires  a  value  given  by  the  Prandtl-Glauert  factor  - - - . 

The  existing  coefficients  for  the  translatory  case  which  are  to  be 
used  in  the  reciprocal  relations  have  been  given  in  various  forms  begin¬ 
ning  with  Possio’s  work.  The  form  in  which  all  the  results  are  commonly 
put  is 

L  =  itpcV2elks  |(ZX  +  iZ2)  (10a) 

and 

M  =  rtpc^e3^8  |(MX  +  iM2)  (lOt) 


where  M  is  the  moment  acting  about  the  quarter- chord  point.  In  order 
to  make  convenient  use  of  the  reciprocal  relations  given  'by  equations  (8a) 
and  (8b),  it  is  first  necessary  to  convert  the  expression  for  lift  given 
by  equation  (10a)  to  a  form  similar  to  the  one  used  in  equation  (9).  When 
the  expressions  (9)  and  (10a)  are  compared,  the  following  correspondence 
is  seen  to  exist  for  compressible  flow: 


*V(k)  = 


_  Z2(k) 


2k 


(11a) 
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and 

Mk)  =  (Ht) 


Sunrraary  of  Available  Flutter  Coefficients  and  Related  Data 
for  Subsonic  Compressible  Flow 

The  two-dimensional  compressible  flutter  coefficients  for  the  real 
and  imaginary  parts  of  the  lift  and  moment  oscillatory  coefficients  for 
sinking  motion  are  given  in  table  I  for  three  Mach  numbers,  M  =  0.5, 

M  =  0.6,  and  M  =  0.7*  This  table  represents  a  summary  of  the  results 
of  three  authors  where  the  flutter  coefficients  by  each  have  been  con¬ 
verted  to  the  Zj_,  Zg,  M-^,  and  M2  form  for  use  in  equations  (10a) 

and  (10b).  The  values  of  Z^_  and  Zg  correspond  to  the  lift  coef¬ 
ficients  indicated  in  equation  (10a),  and  the  values  of  M^  and  M2 

indicated  in  equation  (10b)  correspond  to  the  moment  taken  about  the 
quarter-chord  position.  The  values  of  these  coefficients  were  obtained 
from  three  sources,  and  the  range  of  reduced  frequencies  taken  from  each 
source  is  indicated  in  table  I.  Data  for  large  values  of  reduced  fre-  * 

quency  do  not  appear  to  exist.  The  accuracy  of  these  data  for  reduced 
frequencies  lower  than  approximately  1  is  better  than  the  accuracy  of 
the  data  at  higher  values  of  reduced  frequency  (see  reference  8),  but  a  w 

reasonably  good  solution  for  the  indicial  lift  and  moment  functions  is 
still  possible  since  this  decrease  in  accuracy  at  the  higher  reduced 
frequencies  will  be  shown  to  have  little  effect  on  the  determination  of 
either  the  unsteady  lift  or  unsteady  moment  functions  provided  the  value 
of  the  flutter  coefficients  at  the  infinite  reduced  frequency  is  known. 
Fortunately  these  values  can  be  determined  from  the  initial  value  of  the 
indicial  functions  (see  appendix  A).  If  the  flutter  coefficients  are 
known  up  to  a  reduced  frequency  where  they  have  practically  the  same 
value  as  the  flutter  coefficient  at  the  infinite  reduced  frequency,  then 
they  are  known  for  a  sufficient  range  of  reduced  frequencies  for  deter¬ 
mining  a  reliable  solution  of  the  indicial  functions.  This  condition 
exists  for  M  =  0.7  since  values  of  the  flutter  coefficients  at  a 
reduced  frequency  of  2.5  have  been  found  to  agree  quite  well  with  the 
values  at  the  infinite  reduced  frequency.  However,  values  of  the  flutter 
coefficients  at  M  =  0.5  and  M  =  0.6  are  known  only  for  reduced  fre¬ 
quencies  up  to  1.  Also,  the  flutter  coefficients  at  these  Mach  numbers 
must  be  known  for  higher  values  of  reduced  frequency  than  for  M  =  0.7 
since  the  flutter  coefficients  approach  the  values  at  the  infinite  fre-  * 

quency  at  a  much  slower  rate.  Consequently,  the  only  Mach  number  for 
which  sufficient  data  were  found  available  for  obtaining  the  indicial 
functions  is  M  =  0.7.  * 
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In  order  to  aid  in  the  numerical  solution  of  the  reciprocal  equa¬ 
te  tions,  the  end  points  of  the  Fc(k)  function  should  be  determined.  , 

These  end  points  may  be  determined  independently  of  the  flutter  coef¬ 
ficients.  The  value  of  Fc(0)  which  represents  the  steady-state  lift 

can  be  given  by  the  Prandtl-Glauert  factor  -- — - — When  an  attempt  is 

A  -  m2 

made  to  calculate  this  end  point  from  the  Z^  and  Z 2  coefficients 

given  in  table  I,  an  indeterminate  form  is  found;  however,  the  magnitude 
of  this  indeterminate  quantity  when  properly  evaluated  agrees  very  well 

with  the  factor  The  value  of  Fr(co)  which  corresponds  to  the 

A  -  M2 

value  of  k-^(s)  at  s  =  0 
following  equation 


Numerical  Solution  of  the  Reciprocal  Equation 

*  The  indicial  lift  function  k^(s)  can  be  determined  from  either  of 

equations  (8a)  and  (8b)  where  the  coefficients  Fc(k)  and  Gc(k)  are 

1  determined  from  the  oscillatory  coefficients  by  equations  (lla)  and 

(lib).  It  has  been  noted  that^  of  the  data  presented  in  table  I,  the 
only  Mach  number  for  which  sufficient  data  are  available  for  providing 
a  reasonable  estimate  of  the  indicial  lift  and  moment  on  the  airfoil 
is  M  =  0.7  since  values  of  the  flutter  coefficients  are  not  given  for 
values  of  k  >  1  at  the  other  Mach,  numbers.  A  plot  of  the  Fc(k)  and 

Gc(k)  functions  for  this  Mach  number  is  shown  in  figure  1.  Although  a 
smooth  curve  can  be  drawn  for  the  Fc(k)  function  at  the  higher  reduced 
frequencies,  the  data  for  the  Gc(k)  function  appear  to  be  erratic  in 

this  region.  Consequently,  the  latter  function  is  represented  by  a 
broken  line  at  the  higher  reduced  frequencies.*  Also,  the  numerical 
evaluation  of  the  k-^(s)  function  by  equation  (8b)  is  more  difficult 

than  by  equation  (8a)  because  the  integrand  in  equation  (8b)  is  inde¬ 
terminate  at  k  =  0.  In  order  to  evaluate  this  indeterminate  form,  the 
value  of  the  derivative  of  Gc(k)  at  k  =  0  must  be  known.  Since  the 

function  Gc(k)  is  not  given  in  a  closed  form,  this  derivative  would 

%  be  difficult  to  determine  accurately  in  the  compressible  case.  Because 

of  the  erratic  nature  of  the  Gc(k)  function  at  the  higher  values  of 

.  reduced  frequency  and  because  of  the  presence  of  this  indeterminate  form 


(see  appendix  A)  can  be  determined  by  the 


F  (00)  =  JL 

c  JtM 


(12) 
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of  Gc(k)  at  k  =  0,  equation  (8a)  is  used  throughout  the  ensuing 
analysis  for  determining  the  indicial  functions. 

The  graphical  solution  of  equation  (8a)  may  he  simplified  consider¬ 
ably  by  incorporating  in  the  solution  of  this  equation  the  known  value 
of  the  flutter  coefficient  at  the  infinite  reduced  frequency  by  making 
the  following  substitution  for  the  flutter  coefficient  Fc(k): 


f(k)  =  Fc(k)  -  Fc(~) 


(13a) 


or 


Fc(k)  =  f(k)  +  Fc(co) 


(I3t) 


The  substitution  of  equation  (13b)  into  equation  (lla)  leads  to  an 
alternate  equation  for  k-j_(s): 


kL(s) 


f(k)  sin  ks 
k 


dk 


(14) 


The  graphical  solution  of  k]_(s)  is  thus  simplified  since  the  integral 

in  this  expression  can  be  evaluated  more  readily  than  the  integral  in 
equation  (8a)  because  the  integrand  approaches  zero  much  more  rapidly. 

The  Fc(k)  function  for  M  =  0.7  is  calculated  from  the  flutter 
coefficient  Z2  by  equation  (lla)  and  is  shown  in  figure  1  together 

with  the  value  of  the  asymptote  calculated  by  equation  (12).  The  trans¬ 
formed  function  f(k)  as  obtained  by  equation  (13a)  is  shown  in  fig¬ 
ure  2.  A  plot  of  the  integrand  in  equation  (l4)  for  several  values  of  s 
is  shown  in  figure  3«  These  integrands  were  integrated  by  means  of  a 
planimeter.  The  k-^(s)  function  thus  found  by  means  of  equation  (l4)  is 

shown  in  figure  4  along  with  the  k^(s)  function  for  M  =  0.0.  In  this 

same  figure  a  part  of  an  independent  solution  for  M  =  0.7  given  in 
reference  4  is  also  shown. 

It  will  be  useful  to  fit  the  kj_(s)  function  for  M  =  0.7  given 

in  figure  4  by  some  analytical  function.  Since  the  exponential  function 
has  a  simple  operational  equivalent  and  has  been  found  convenient  in 
approximating  the  k-j_(s)  function  at  M  =  0.0  (see  reference  ll),  a 
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limited  series  of  such  functions  was  chosen  to  approximate  the 

k^(s)  function  at  M  =  0.7*  The  function  found  to  fit  this  curve  quite 

well  is 


( s )  =  1.4(l  -  0.364e"0, °536s  _  o.405e"°‘  +  0.4l9e"°*9°2s)  (15) 


and  is  also  shown  in  figure  4.  The  corresponding  approximate  expressions 
for  the  harmonically  oscillating  airfoil  can  he  found  from  equation  (15) 
in  a  similar  manner  shown  in  reference  11.  The  expressions  for  Fc(k) 
and  Gc(k)  are 


F_(k)  =  1.4 


0.364k' 


(0.0536) 2  +  k2 


0. 405k‘ 


(0.357) 2  +  k2 


0. 4lQk‘ 


2  2 
(0.902)  +  k 


(16a) 


and 


Gc(k)  =1.4 


Both  of  these  expressions  are  plotted  in  figure  1  for  comparison  with 
the  data  obtained  from  the  flutter  coefficients.  Note  that  the  approxi¬ 
mate  expression  for  the  Fc(k)  function  is  in  fairly  good  agreement  for 
all  of  the  k  values  considered.  However,  relatively  poor  agreement 
exists  for  the  Gc(k)  function  at  the  higher  values  of  k.  This  com¬ 
parison  is  to  be  expected  since,  as  was  previously  cited,  the  Gc(k)  func¬ 
tion  obtained  from  the  flutter  coefficients  appears  to  be  quite  erratic 
at  the  higher  values  of  k. 

The  solution  of  the  reciprocal  equation  has  thus  far  been  shown  for 
the  calculation  of  the  k-^(s)  function.  Since  the  reciprocal  equations 

are  also  applicable  to  the  determination  of  the  unsteady  moment  rn^(s) 
due  to  a  sudden  change  in  airfoil  vertical  velocity,  the  rn^(s)  function 
taken  about  the  quarter- chord  point  can  be  shown  to  be  given  by 


0.01951k 


0.l446k 


(0.0536) 2  +  k2 


2  2 
(0.357)  +  k 


0.3779k 


2  2 
(0.902)  +  k 


(16b) 
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(s)  .  £  r  MW_°m  fa  a*  (17) 

^  « J  0  k 


where  M(k)  is  considered  herein  to  he  determined  from  the  flutter  coef¬ 
ficient  Mg  hy  the  following  equation: 

M(k)  -  (18) 


A  plot  of  the  function  M(k)  for  M  =  0. 7  is  shown  in  figure  5  together 
with  the  value  of  M(“)  determined  by 


M(°°) 


1 

2nM 


(19) 


In  a  manner  similar  to  that  shown  for  the  unsteady- lift  case,  the  func¬ 
tion  M(k)  is  then  transformed  in  order  that  equation  (17)  can  be 
evaluated  graphically  by  the  following  substitution: 


m(k)  =  M(k)  -  M(») 


(20a) 

* 


or 


M(k)  =  m(k)  +  M(°°) 


(20b) 


A  plot  of  the  function  m(k)  for  M  =  0.7  is  shown  in  figure  6.  The 
substitution  of  equation  (20b)  into  equation  (17)  leads  to  the  expression 

m1(s)  =  M(oo)  +  £  f  m(k)  sxn  ks  ^  (21) 

“JO  k 

Plots  of  the  integrand  in  equation  (21)  are  shown  in  figure  7  for  several 
values  of  s,  and  the  m^(s)  function  evaluated  by  this  equation  is  shown 

in  figure  8. 

If  the  indicial  lift  function  k^s)  and  the  indicial  moment  func¬ 
tion  m-^(s)  are  known,  the  variation  of  indicial  center-of-pressure 
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^  location  in  percent  chord  can  be  calculated  by  means  of  the  following 

equation: 


xcp(s) 


=  25  ml^s^100 

kx(s) 


-(22) 


A  plot  of  the  center- of-pres sure  location  is  given  in  figure  9  and.  for 
■comparison,  the  case  for  M  =  0.0  is  also  shown. 


Approximate  Considerations  for  the  Indicial  Lift  Due  to 
Penetration  of  a  Sharp-Edge  Gust 

Because  of  the  limited  number  of  flutter  coefficients  available  for 
flap  motions,  as  noted  in  the  introduction,  the  determination  of  the 
indicial  lift  function  for  penetration  k2(s)  by  use  of  the  reciprocal 

equations  is  not  possible  at  the  present  time.  An  alternate  and  perhaps 
semirational  method  for  determining  a  part  of  the  penetration  func¬ 
tion  k2(s)  from  the  k^s)  function  on  the  basis  of  the  relations 
between  these  two  functions  for  the  incompressible  case  is  suggested  as 
follows.  The  relations  between  the  k-^(s)  and  k2(s)  functions  for  the 

incompressible  case  are  given  in  reference  3  and  are  as  follows: 


k2(s)  ~'*fQ  kl(S  "  Sl)|/2  _1Bi  dsl  +  ilM2  "  s)  (0  <  s  <  2)  (23a) 


k2(s)=ij^  kl(S  "  Sl)^F~  dsl  (s  >  2)  (23b) 

The  first  term  in  equation  (23a)  is  associated  with  only  the  circulatory 
part  of  the  k-|_(s)  function  while  the  second  term  of  the  equation  arises 

from  the  impulsive  term  in  the  k-^s)  function,  and  hence  is  associated 

with  the  apparent-mass  effects.  This  second  term  is  known  not  to  apply 
*r  for  compressible  flow,  since  for  the  kx(s)  function  the  apparent-mass 

effect  is  no  longer  impulsive  in  character  but  rather  is  time  dependent. 
Therefore  this  equation  cannot  be  applied  to  the  compressible  case.  With 
regard  to  the  applicability  of  equation  (23b),  it  is  assumed  that  this. 
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equation  can  evaluate  the  circulatory  part  of  the  lift  for  the  k2(s)  func¬ 
tion  provided  that  only  the  circulatory  part  of  the  k^(s)  function  is 
used.  Although  the  k-^s)  function  calculated  herein  for  compressible 

flow  contains  hoth  the  circulatory  and  noncirculatory  lifts,  the  following 
argument  may  he  used  as  a  basis  for  the  assumption  that  for  values  of  s 
greater  than  approximately  4,  the  k-^(s)  function  can  he  considered,  for 

practical  purposes,  to  he  due  to  circulation  only.  If  the  assumption  is 
made  that  the  lift  due  to  apparent-mass  effects  and  the  lift  due  to  circu¬ 
lation  act  always  at  i-- chord  position  and  chord  position,  respectively, 

the  moment  function  m1(s)  shown  in  figure  8  must  arise  from  the  apparent- 

mass  effects  only.  This  figure  shows  that  for  practical  purposes  the 
moment  function,  referred  to  the  quarter-chord  position,  has  decayed  to  a 
negligible  value  for  s  greater  than  approximately  4.  Thus  the  portion 
of  the  k-]_(s)  function  in  figure  4  for  s  greater  than  approximately  4 

may  he  associated  with  the  lift  due  to  circulation  only.  Inspection  of 
equation  (23b)  shows  that,  if  the  evaluation  of  the  k2(s)  function  is 
confined  to  values  of  s  greater  than  6,  only  the  portion  of  the 
k-|_(s)  function  for  s  greater  than  approximately  4  (the  portion  attributed 

to  circulation)  will  he  used.  A  numerical  solution  of  equations  (23a) 
and  (23b)  is  given  in  appendix  B.  The  k2(s)  function  for  M  =  0.0  was 

calculated  by  using  the  kx(s)  function  for  M  =  0.0  shown  in  figure  4 

as  a  means  for  estimating  the  accuracy  of  a  numerical  solution  given  by 
equation  (B9)  in  appendix  B.  The  values  obtained  were  within  1.5  percent 
of  the  known  Kiissner  function.  For  the  compressible  case  the  k2(s)  func¬ 
tion  for  s  greater  than  approximately  6  was  evaluated  by  utilizing  the 
k-j^s)  function  for  M  =  0. 7.  The  results  of  this  calculation  are  shown  in 

figure  10.  Although  the  solution  for  the  compressible  case  is  approximate, 
it  is  interesting  to  note  that  for  s  >6,  where  the  solution  was  reasoned 
to  be  applicable,  the  growth  of  lift  is  less  rapid  for  the  con^ressible  case 
than  for  the  incompressible  case.  This  phenomenon  would  be  expected  since 
the  growth  of  lift  due  to  the  kx(s)  function,  which  was  used  to  determine 

the  k2(s)  function,  shows  this  same  effect  as  evidenced  by  figure  4. 


POSSIBLE  EFFECTS  ON  GUST  LOAD  FACTOR 


The  fact  that  the  growth  of  lift  for  the  k-^s)  and  k£(s)  functions 

is  less  for  the  M  =  0.7  case  than  for  the  M  =  0.0  case  may  have  some 
significance  in  the  determination  of  the  response  of  an  airplane  to  gusts. 

The  current  design  procedure  for  the  loads  due  to  gusts  on  an  airplane 
utilizes  the  two-dimensional  incompressible- flow  k^(s)  and  k2(s)  functions 
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and  an  over-all  approximate  correction  for  compressibility.  If  in  the 
calculation  of  the  gust  load  factor  at  about  M  =  0.7  the  corresponding 
k-j_(s)  and  ^(s)  functions  for  M  =  0.7  were  used  instead  of  the 

incompressible  k-^(s)  and  k2(s)  functions  with  the  over-all  approxi¬ 
mate  correction  for  compressibility,  then  the  resulting  load  factor 
would  probably  be  less  than  that  found  by  using  the  incompressible 
kj_(s)  functions.  This  further  suggests  the  desirability  of  obtaining 

the  k2(s)  function  for  the  complete  range  of  s  and  of  obtaining  the 

k2(s)  functions  for  other  values  of  Mach  number  so  that  a  greater 

insight  may  be  had  as  to  the  effect  of  compressibility  on  the  gust  load 
factor. 


CONCLUDING  KEMAEKS 


The  reciprocal  equations  which  relate  the  incompressible  circulatory 
lift  to  the  lift  due  to  harmonic  oscillations  have  been  modified  to 
include  the  noncirculatory  lift  due  to  apparent-mass  effects.  While  the 
lift  due  to  apparent-mass  effects  for  the  indicial  case  is  known  to  be 
impulsive  for  the  incompressible  case,  the  lift  due  to  apparent-mass 
effects  for  the  corresponding  compressible  case  was  found  to  be  finite 
and  a  time-dependent  function. 

The  reciprocal  equations  are  evaluated  numerically  for  determining 
the  indicial  lift  and  moment  on  an  airfoil  experiencing  a  sudden  acqui¬ 
sition  of  vertical  velocity  at  M  =  0.7*  Lack  of  sufficient  flutter 
coefficients  prevents  the  calculation  of  these  functions  at  other  Mach 
numbers. 

The  growth  of  lift  to  the  steady- state  value  on  an  airfoil  experi¬ 
encing  a  sudden  acquisition  of  vertical  velocity  was  found  to  be  less 
rapid  for  the  compressible  case  than  for  the  incompressible  case. 

An  approximate  calculation  is  made  for  a  portion  of  the  indicial- 
lift  function  due  to  penetration  of  a  sharp-edge  gust  at  M  =  0.7*  It 
is  noted  that  the  growth  of  lift  to  the  steady- state  value  for  this 
indicial  lift  function  was  also  less  rapid  than  for  the  incompressible 
case.  As  a  consequence  of  these  phenomena,  the  calculation  of  the  gust 
load  factor  at  high-subsonic  Mach  numbers  utilizing  the  two-dimensional 
incompressible  indicial  lift  functions  and  an  over-all  correction  for 
compressibility  such  as  the  Prandtl-Glauert  factor  might  be  conservative. 

Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  July  30,  1951 
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APPENDIX  A 

« 


PROOF  FOR  lim  k-.  (s)  =  F,,(oo) 
s — >  +0 

The  function  k-^(s)  can  he  expressed  as  a  function  of  Fc(k)  by 
the  following  equation: 


_  2  r  Fc(k>  »i.  k3 

1  n  J  0  k 


dk 


(Al) 


provided  Fc(k)  is  a  bounded  and  well-behaved  function.  The  value  of 

k^(s)  as  s  tends  to  zero  positively  may  by  evaluated  by  taking  the 
limit  of  both  sides  of  equation  (Al)  as  s — >+0: 


lim  kj ( s ) 
s  -^>+0 


lim 
s — >+0 


2  r  Fc(k)  sin 

*J0 


dk 


(A2) 


Before  the  limits  in  equation  (A2)  are  evaluated,  it  is  convenient  to 
make  the  following  substitution  for  Fc(k): 


Fc(k)  =  f(k)  +  FCH 


(A3) 


Substitution  of  this  expression  for  Fc(k)  into  equation  (A2)  leads  to 
the  following  expression  for  k-^(s): 


lim  k^(s) 
s  — >  0 


(A4) 
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Although  the  second  integral  in  equation  (A4)  can  he  evaluated  readily 
*  to  give 


lim  P  ■  dk  =  ~  (A5) 

s — >+0^0  k  2 

the  evaluation  of  the  first  integral  in  equation  (a4)  is  not  readily 
evident;  however  it  can  he  shown  (see,  lemma)  that 


lim  r*W.  an-.  0 

s — )0 0  k 


(A  6) 


since 


lim  f(k)  =  0 

k — ^  00 


As  a  consequence  of  equations  (A5)  and  (A6)  the  value  of  the  limit  in 
equation  (A4)  is 


lim  k-Js)  =  Fc(»)  (A?) 

s — >+0 


Lemma: 

Prove  that 


lim 
s — >+0 


f(k)  sin  ks 


dk  =  0 


(A8) 


provided  f(k)  is  a  hounded  and  integrahle  function  within  the  limits  of 
integration  and 


lim  f (k)  -  0 

k — 00 
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If  the  following  substitution  is  made  for  ks 

« 


ks  =  x 


an  alternate  expression  for  equation  (A8)  is  had 


lim  f  f(|)  dx  =  0  (A9) 

s — )+0  ^0 


Since  f(f-)  is  hounded  and  integrable  within  the  limits  of  integration 
and  since  the  integral 


is  absolutely  convergent,  then  equation  (A9)  is  uniformly  convergent 
for  s  =  0  within  the  limits  of  integration.  Consequently,  there  is 
a  positive  number  X  such  that 


^sin  x  dx 

r  f(\ 

Jx  U 

x 

» 


(A10) 


where  6  is  an  arbitrary  number  chosen  as  small  as  desired  and  the 
number  X  is  independent  of  the  value  of  s.  Also  a  number  xQ,  inde¬ 
pendent  of  s,  may  be  chosen  so  that 


S0  0  f©- 


sm  x 


dx 


6 


(All) 


However, 


/ x  *©*4* 


^2^  dx  +  f 
X 


(5)/. 


sm  x 


dx 


x 
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vhere 


xQ<J<X 


It  can  be  shown  that 


therefore 


<  it 


(q.  >  p  >  0) 


<  2*f(£) 


thus  a  value  of  S  can  he  chosen  as  small  as  desired  so  that 

X 
x 


/  <D- 


sm  x 


dx 


<*• 


(s  %  S)  (A12) 


Combining  the  results  given  by  equations  (A10),  (All),  and  (A12)  yields 
the  equation 


r“ 

/  f(^ 

c\sin  x  dx 

Jo  . 

W  X 

<h+h+^s 

3  3  3 


(s  ^  S)  (A13a) 


thus,  in  the  limit  as  S — >0  positively 


lim 
s — >+0 


f/x\sin  x 
Vs/  x 


dx  =  0 


(A13b) 


or 


lim 
s — >  +0 


f(k)sin  ks  dk 
k 


(A13c) 
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APPENDIX  B 


NUMERICAL  EVALUATION  OF  INCOMPRESSIBLE-FLOW  INDICIAL  LIFT  FUNCTION 
DUE  TO  PENETRATION  OF  A  SHARP-EDGE  GUST 


The  indicial  lift  function  due  to  penetration  of  a  sharp-edge 
gust  kgCs)  expressed  in  terms  of  the  indicial  lift  function  k^s) 

is  given  by  equations  (23a)  and  (23b).  Although  a  numerical  evaluation 
of  the  first  term  in  equation  (23a)  may  be  performed  with  the  aid  of 
Simpson's  rule,  the  numerical  evaluation  of  the  integral  in  equation  (23b) 
is  difficult  in  that  the  integrand  is  infinite  at  s ^  =  2.  In  order  to 

overcome  this  difficulty,  a  parabolic-arc  segment  can  be  fitted  to  the 

function  k-^s  -  s^,  between  the  limits  s-j_  =  ^  an(*  si  =  2,  Because 

of  the  similarity  between  equations  (23a)  and  (23b),  the  numerical  solu¬ 
tion  of  both  equations  can  be  written  in  one  equation  with  the  aid  of 
matrix  notation.  For  simplicity,  the  numerical  solution  is  calculated 
at  integral  values  of  s.  However,  for  increased  accuracy  the  interval 
of  integration  is  taken  at  quarter  multiples  of  s.  In  view  of  these 
considerations,  the  numerical  solution  for  the  kp^)  function  may  be 

performed  in  three  steps: 

(l)  Integration  of  integrand  for  the  limits  s-^  =  0  to  s^  =  1 

With  the  aid  of  Simpson* s  integrating  factors  the  following  integral 
is  evaluated,  a  quarter  multiple  of  s  being  assumed  for  the  interval 
of  integration: 


*  f0  ki(e  - 

(2)  Integration  of  integrand  for 


the  limits 


=  1  to 


(Bl) 


In  a  manner  similar  to  that  shown  for  the  evaluation  of  the  integral 
in  step  1,  Simpson1 s  integrating  factors  can  be  utilized  to  calculate  the 

same  integral  over  the  limits  3-^  =  1  to  s-j_  =  ^  as  follows: 


* 


# 
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1 — 1 
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1 

C\J 

1  14 
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3  ki(s  -  ?) + 1  "  i 


(3)  Integration  of  integrand  for  the  limits  si  =  2  si  =  2 

As  noted  previously  the  integration  for  these  limits  is  not  possible 
by  means  of  the  simple  Simpson’s  factors  since  the ■ integrand  is  infinite 
at  sq_  =  2.  With  the  change  in  variable  s-^  =  2  -  a,  a  parabolic-arc 

segment  is  fitted  to  the  function  k-^s  -  s.jJ  between  the  limits  s^  - 
to  s^  =  2  by  use  of  the  following  equation: 


k1(s~2  +  <j)*'A  +  Bcr  +  Ccr2  (B3) 

A  set  of  simultaneous  equations  relating  the  function  k-^s  -  2  +  cr)  with 

the  constants  A,  B,  and  C  in  the  parabolic-arc  segment  may  be  obtained 

by  evaluating  equation  (B3)  at  the  points  a  =  0,  cr  =  i,  and  a  =  1  as 

4  2 

follows:  at  a  =  0, 


kx(s  -  2)  =  A 


at  cr  =  i 


k-.  (b  -  — ^  =  A  +  i  B  +  —  C 

H  k)  4  16 


at  a  =  ± 


;i(s  - 1)  ■ A  +  5 B  + 1 0 


OOlCVJ 
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When  equations  (b4) ,  (B5),  and  (B6)  are  solved  simultaneously  for 
the  constants  A,  B,  and  C,  the  resulting  expressions  are 


A  =  k-Js  -  2) 


B  =  -6^(3  -  2)  +  Sk-^s 
C  =  8k1(s  -  2)  -  l6kx (a 


5)  -  - §) 

£W-f) 


(B7) 


* 


* 


If  the  expression  for  k-j_(s  -  2  +  a)  given  by  equation  B(3)  is  substi¬ 
tuted  into  the  integrand  under  consideration  together  with  the  values 
of  A,  B,  and  C  determined  by  equations  (B7),  the  integral  can  be 
evaluated  in  closed  form  in  terms  of  the  unsteady-lift  function  k^(s  -  s-]J 

from  s^  =  ~  to  s-j_  =  2.  The  results  of  this  closed- form  integration  are 
as  follows: 


In  figure  10  the  solutions  for  the  k^s)  function  are  given 

for  M  =  0.0  and  M  =  0.7.  No  discernahle  difference  was  present  "when 
the  calculated  values  for  M  =  0. 0  were  compared  with  the  Kussner 
function. 
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IMACA 


functions  due  to  translator^-  oscillation 


Figure  3.—  Plots  of  the  integrand  in  equation  (lit.)  for  determining  the 
.  indicial  lift  k-i(s)  at  M  =  0.7. 
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